Tensor network states for quantum spin ladders 



O 
(N 



Sheng-Hao Li,' Yao-Heng Su,' Yan-Wei Dai,' and Huan-Qiang Zhou' 

' Centre for Modern Physics and Department of Physics, 
Chongqing University, Chongqing 400044, The People's Republic of China 

We have developed an efficient tensor network algorithm for spin ladders, which generates ground-state wave 
functions for infinite-size quantum spin ladders. The algorithm is able to efficiently compute the ground-state 
fidelity per lattice site, a universal phase transition marker, thus offering a powerful tool to unveil quantum 
many-body physics underlying spin ladders. To illustrate our scheme, we consider the two-leg and three-leg 
Heisenberg spin ladders with staggering dimerization. The ground-state phase diagram thus yielded is reliable, 
compared with the previous studies based on the density matrix renormalization group. Our results indicate that 
the ground-state fidelity per lattice site successfully captures quantum criticalities in spin ladders. 
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Introduction. Tensor networks (TN) provide a convenient 
means to represent quantum wave functions in classical sim- 
ulations of quantum many-body lattice systems, such as the 
matrix product states (MPS) Uj-lsl] in one spatial dimension 
and the projected entangled-pair state (PEPS) [6,-8] in two 
and higher spatial dimensions. The development of various 
numerical algorithms in the context of the TN representations 
has led to significant advances in our understanding of quan- 
tum many-body lattice systems in both one and two spatial 
dimensions i3l 4l7ll . Lying between quantum lattice systems 
in one and two spatial dimensions, spin ladders have attracted 
a lot of attention, due to their intriguing critical properties. 
Given the importance of spin ladder systems in condensed 
matter physics, it is somewhat surprising that no efforts have 
been made to develop any efficient algorithm in the context of 
the TN representations. 

This paper aims to fill in this gap. The algorithm gener- 
ates efficiently ground-state wave functions for infinite-size 
spin ladders. In addition, it allows to efficiently compute the 
ground-state fidelity per lattice site, a universal phase transi- 
tion marker, thus offering a powerful tool to unveil quantum 
many-body physics underlying spin ladders. In fact, as ar- 
gued in Refs. [18-23], the ground-state fidelity per lattice site 
is able to capture drastic changes of the ground-state wave 
functions around a critical point. To illustrate our scheme, 
we consider the two-leg and three-leg Heisenberg spin lad- 
ders with staggering dimerization. The ground-state phase di- 
agram thus yielded is reliable, compared with the previous 
studies a24l 12511 based on the density matrix renormalization 
group (DMRG) | |26ll . Our results indicate that the ground-state 
fidelity per lattice site successfully captures quantum critical- 
ities in spin ladders. 

Tensor network representation for spin ladders. Let us de- 
scribe the TN representation suitable to describe a ground- 
state wave function for an infinite-size spin ladder Suppose 
the Hamiltonian is translationally invariant under shifts by ei- 
ther one or two lattice sites along the legs: H - Yi{i,a) h(i,a), 
with the (/, Q')-th plaquette Hamiltonian density h(i^a) acting 
on sites / and {i + 1) along the a-th and {a + l)-th legs. 
Here, (/, a) runs over all the possible plaquettes by taking 

e {-oo, ■ ■ ■ , +oo), and a - l,---,n-l, with n being the 



number of the legs. Assume that the TN representation for a 
wave function enjoys the translational invariance under shifts 
by two lattice sites along the legs. In the following, we focus 
on a detailed description for a two-leg spin ladder, with a brief 
discussion for an n-leg ladder system. 

For an infinite-size two-leg spin ladder system, we need 



only four different four-index tensors A^^^, B^^^ 
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to Store the wave function. Here, A 
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made of complex numbers labeled by one physical index s 
and four inner bond indices {, r, u and d, where i = 1, d, 
with d being the dimension of the local Hilbert space, and 
{!, r, u, d - 1,...,D, with D being the bond dimension. A 
four-index tensor A^^^ is visualized in FiglTJi), with a similar 
pictorial representation for the tensors B^^^, O^^^, and A 
TN representation for the ground-state wave function is shown 
for an infinite-size two-leg spin ladder in FiglTJii). There are 
two different but equivalent choices of the unit cell for such 
an infinite-size TN: one is chosen as A, B, D, and C clockwise 
if / is even, the other is chosen as B, A, C, and D clockwise if 
/ is odd, see FiglTJiii)- 

Now let us turn to the computation of the norm for a quan- 
tum state wave function. To this end, we introduce double 
tensors a^f^, b(fj, q^;}, and df^^, with £ = ((,('), f = (r,r'), 
u = (u, u ), and d = (d, d ). They form from the four-index 
tensors Aj.^^, B^^^, C'J^^^, and DJ^^^, and their complex conju- 
gates, see Figdiiv) for a pictorial representation of the double 
tensors. With these double tensors, the TN for the norm of a 
wave function is shown in FiglTJv). Again, we have two dif- 
ferent but equivalent choices for the unit cell of the norm TN: 
one is a, b, d, and c clockwise if / is even, the other is b, a, c, 
and d clockwise if / is odd, see Figdivi). 

The expectation value of an operator acting on a plaque- 
tte, such as the ground-state energy per unit cell, also admits a 
TN representation, which absorbs the operator acting on a pla- 
quette for an infinite-size spin ladder system. For a randomly 
chosen initial state |i/^o)j the energy is expressed as. 
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For different choices of the unit cell, we get two different but 
equivalent forms of the zero-dimensional transfer matrix E 
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FIG. 1 : (Color online) (i) Four-index tensor A-'.^ used to represent a 
TN representation for the ground-state wave function for an infinite- 
size two-leg spin ladder, with s being a physical index, /, r, and d 
denoting the inner indices, (ii) The pictorial representation for a TN 
state |i/r> with leg and rung bonds, which are used to absorb an opera- 
tor acting on the ;-th plaquette. (iii) Two different choices of the unit 
cell for an infinite TN state, made of four four-index tensors A, B, C, 
and D. (iv) A double tensor a-ffj is formed from the four-index tensor 
A^^^ and its complex conjugate (A*)^, , with E = ), f = (r, r ), 

and d = (d,d ). (v) The TN representation for the norm of a ground- 
state wave function in an infinite-size spin ladder, (vi) Two different 
choices of the unit cells for the norm tensor net work, consisting of 
four double tensors a, b, c, and d. 




FIG. 3: (Color online) The contribution to the energy gradient for an 
infinite-size spin ladder consists of three parts: (i) the hole cell with 
the tensor A^^^ removed and the Hamiltonian cell with the Hamilto- 
nian density acting on a plaquette locate on the same cell; (ii) the 
hole cell locates on the right hand side of the Hamiltonian cell; (iii) 
the hole cell locates on the left hand side of the Hamiltonian cell. In 
the latter two cases, there are m cells between the hole and Hamilto- 
nian cells, where m e (0, 1, 2, 3, ■ ■ • ). Here, the hole cell is visualized 
in (iv), with the tensor A^^.^ removed. 
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FIG. 2: (Color online) The ground-state energy per unit cell is com- 
puted for translation-invariant spin ladders, (i) The transfer matrix E 
for an infinite-size norm tensor network, which is constructed from 
four double tensors a^fj, bff^, Cffn, and dff^,, with { = (l,E ),f = (r, r ), 
u = (u,u ), and d = (d,d ). (ii) The dominant left and right eigen- 
vectors Vi and Vr of the transfer matrix E. (iii) A unit cell with the 
Hamiltonian density hAsoc acted on the plaquette. (iv) The ground- 
state energy per unit cell is computed from the eigenvectors V/., Vr, 
four four-index tensors A! ,, Bl ,, Ci , and Dl . 

constructed from four double tensors a-^fj, b~fjj, djj.j^, and Cf,,j,, 
one of them is shown in Fig.|2fi). The dominant left and right 
eigenvectors of the transfer matrix E constitute the environ- 
ment tensors, visualized in Fig. Etii). This enables us to ab- 
sorb an operator acting on the ;-th plaquette A, B, D, and C 
clockwise, if / is even, as shown in Fig. IJiii), and compute 
the energy per unit cell, as shown in Fig.|2liv). The same pro- 
cedure may be used to compute the energy per unit cell for an 
operator acting on the i-th plaquette B, A, C, and D clockwise, 
if / is odd. 

To update the TN representation, we compute the energy 
gradient with respect to four-index tensors: 

dE diUHmidAl^ dm<Pa)ldAl^ 

= E ■ . (2) 

5A^^^ <<Ao|iAo) (lAoliAo) 

Here, a four-index tensor A^^^ is used to explain how to effi- 
ciently evaluate the energy gradient in the context of the tensor 
network representation for an infinite-size two-leg spin ladder. 



with the details visualized in Fig. [3] Notice that the contribu- 
tions to the energy gradient come from three parts: (i) the hole 
cell with the four-index tensor A^^^ absent and the Hamilto- 
nian cell with the Hamiltonian density sandwiched locate on 
the same cell; (ii) the hole cell locates on the right hand side 
of the Hamiltonian cell; (iii) the hole cell locates on the left 
hand side of the Hamiltonian cell. In both cases (ii) and (iii), 
there are m cells between the hole cell and the Hamiltonian 
cell, where m € (0, 1, 2, 3, ■ • •). As such, the four-index tensor 
A^^^j is updated as follows, 

BE 

(3) 

where S denotes the step size during updating. We stress that, 
for a two-leg spin ladder, we should update four different four- 
index tensors A^^.^^, B^^^, C^^.^^, and D^^^^ simultaneously. 

The above updating procedure yields new tensors A^^^, B^^.^, 
C^i^.^, and D^^^ for a two-leg spin ladder Repeating this pro- 
cedure until the ground-state energy per unit cell converges, 
we anticipate that the system's ground-state wave function is 
generated in the TN representation. 

For a three-leg spin ladder, one should introduce four dif- 
ferent four-index tensers Ai ,, Bi ,, Ei , and Fl , and two 
different five-index tensers C^^„^ and £'^„,j. Similarly, more 
tensors are needed for a multi-leg spin ladder However, the 
algorithm is applicable to a multi-leg spin ladder, as long as 
the memory is sufficient to store the TN tensors. 

The models. As an illustration, we test the algorithm on the 
infinite-size two-leg and three-leg Heisenberg ladder systems 
with staggering dimerization. 

The two-leg and three-leg Heisenberg spin ladders are, re- 
spectively, described by the Hamiltonian 

+ J±Y,^u-Sia, (4) 

11=1,2 ! i 
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TABLE I: The extrapolated infinite-size ground-state energy per site 
from finite-size spin- 1/2 ladders in Ref. |27] vs. our ground-state 
energy per site for the infinite-size two-leg Heisenberg ladder. 



and 

H- ^ ''^Ji,aSi,a-Si+\,a + J±''^{Si^l-Si^2 + SixSi,i) (5) 
a= 1,2,3 ! i 

where denotes the spin- 1/2 Pauli operator at site / on the 
Qf-th leg, and J^^ is the exchange interaction coupling along the 
rungs. 

In order to test our algorithm, we first consider a two-leg 
ladder Choose the exchange interaction coupling constant 
Ji_ e 0.2,0.4,0.6,0.8, 1.0, and the coupling constant in each 
chain - J ( J - 1) if / -H or is odd, - J' if i + a is even, 
with r e 0, 0.2, 0.4, 0.6, 0.8, 1.0. In Table U we list our simu- 
lation results for the ground-state energy per site, for different 
values of and J' , with the truncation dimension up to 6, 
against the extrapolated infinite-size ground-state energy per 
site from finite-size spin- 1/2 ladders in Ref.|27]. The fact that 
they matches very well demonstrates that our TN algorithm 
for spin ladders is reliable. 



Second, we focus on critical points of the ladders with the 
staggered dimeiization 7,> = J{\ + (-l)'^"^], which are the 
exchange interaction couplings along the a-th leg for the two- 
leg ladder {a - 1,2) and the three-leg ladder {a - 1,2, 3). In 
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FIG. 4: (color online) The ground-state fidelity per lattice site 
d(J^i, J j^2), as a function of Jj^i and Jj^2 for the two-leg Heisen- 
berg ladder with staggering dimerization. Upper panel: A two- 
dimensional fidelity surface embedded in a three-dimensional Eu- 
clidean space. A continuous phase transition point J±c - 1.24 is 
identified as a pinch point (J±c,J±c) on the fidelity surface, as ar- 
gued in Refs. Il8l42in. Here, we have taken the truncation dimension 
D = 6. Lower panel: The contour plot of the fidelity per lattice 
site d(J^[, J j^2), on the (7j^i , yx2)-plane, for the two-leg Heisenberg 
ladder with staggering dimerization. 



addition, we choose 6 = 0.5, and the coupling constant J to 
be unity (J = 1). To this end, we need to compute the fidelity 
per lattice site. 

The ground-state fidelity per lattice site. As argued in 
Refs. 1118142 in . the ground-state fidelity per lattice site is a uni- 
versal marker to detect a quantum phase transition: a phase 
transition point is characterized by a pinch point on the fidelity 
suiface. 

Consider the Heisenberg ladders with staggering dimeriza- 
tion. We choose the exchange interaction coupling along the 
rungs as a control parameter For two different ground 
states, |i/'(7j_i)) and |i/'(7j_2)) coiTesponding to two different 
values and 7^2 of the control parameter 7^, the ground- 
state fidehty F(7j.i,7_l2) - Ki/'(-^-l2)I'/'(-^-li))I asymptotically 
scales as F(7j.i,7j.2) ~ d(J j_i , J j_2)'^ , with the system 
size. Here, d(J±\,J±2) is the scaling parameter, introduced 
in Refs. flI8l - 2011 for one-dimensional quantum lattice systems 
and in Ref. [2 1] for two and higher-dimensional quantum lat- 
tice systems; it characterizes how fast the fidelity between two 
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FIG. 5: (color online) The ground-state fidelity per lattice site 
d{J^\ , J±2), as a function of J^i and 7^2 for the three-leg Heisen- 
berg ladder with staggering dimerization. Upper panel: A two- 
dimensional fidelity surface embedded in a three-dimensional Eu- 
clidean space. A continuous phase transition point J±c - 0.96 is 
identified as a pinch point (J±c,J±c) on the fidelity surface, as ar- 
gued in Refs. |18l42ll]. Here, we have taken the truncation dimension 
D = 6. Lower panel: The contour plot of the fidelity per lattice 
site d{J ^i,J ^2), on the (7x1. 7x2)-plane, for the two-leg Heisenberg 
ladder with staggering dimerization. 



panel, a contour plot is shown for the fidelity per lattice site 
on the (7x1, 7x2)-plane. We stress that no significant shifts are 
observed for the pinch point, when the truncation dimension 
is increased up to 6. Therefore, we conclude that a continu- 
ous phase transition takes place at 7xc - 1 -24, which is very 
close to earlier results from the mean-field theory ll28ll , exact 
diagonalization |29] and DMRG (H. 

Similar to the two-leg Heisenberg spin- 1/2 ladder, we 
plot a two dimensional fidelity surface embedded in a three- 
dimensional Euclidean space, namely, the ground-state fi- 
delity per lattice site d(Ju,J±2) as a function of 7xi and 7j.2, 
for the three-leg Heisenberg spin- 1/2 ladder with staggering 
dimerization in Fig.|5] It yields reliable results, with only the 
truncation dimension up to 6. A continuous phase transition 
point 7xf - 0.96 is identified as a pinch point (7j.c, J±c) on the 
fidelity suiface, consistent with the previous results from the 
DMRG method (Hill. 

Conclusions. We have developed an efficient TN algorithm 
to compute ground-state wave functions for infinite-size quan- 
tum spin ladders. Our investigation lends fuither support to 
the observation that the ground-state fidelity per lattice site 
is able to characterize critical phenomena in quantum many- 
body systems. It also demonstrates that the developed TN 
algorithm for spin ladders is efficient to compute the fidelity 
per lattice site. 
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ground states goes to zero when the thermodynamic limit is 
approached. Physically, the scaling parameter d(Ju,J±2) is 
the averaged fidelity per lattice site, 

lnd(J_Li,J±2)= lim , (6) 

which is seen to be well defined in the thermodynamic limit. It 
satisfies the properties inherited from the fidelity /^(7j^i, 7^2): 
(i) normalization (i(7x,/±) = 1; (ii) symmetry d(Ju,J±2) - 
d(J±2, and (iii) range < d{J±i,J±2) ^ 1. 

We emphasize that the TN representation of the system's 
wave functions generated from the algorithm makes it efficient 
to compute the ground-state fidelity per lattice site for spin 
ladders. 

In FiglH we plot the ground-state fidelity per lattice site 
for the two-leg Heisenberg spin- 1/2 ladder with staggering 
dimerization. A two-dimensional fidelity suiface embedded 
in a three-dimensional Euclidean space is shown in the upper 
panel, with a pinch point at (1.24,1.24), implying that a con- 
tinuous phase transition occurs at 7xc - 1.24. In the lower 
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